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ABSTRACT 

We have determined the central velocity dispersions and surface brightness profiles 
for a sample of powerful radio galaxies in the redshift range 0.06 < z < 0.31, which 
were selected on the basis of their young radio source. The optical hosts follow the 
fundamental plane of elliptical galaxies, showing that young radio sources reside in 
normal ellipticals, as do other types of radio galaxies into which these objects are 
believed to evolve. As young radio sources are relatively straightforward to select and 
the contributions of the AGN light to the optical spectra are minimal, these objects can 
readily be used to study the evolution of the fundamental plane of elliptical galaxies 
out to z=l, independently of optical selection effects. 

The black hole masses, Afh^,, of the objects in our sample have been determined 
using the tight empirical relation of Mi^h with central velocity dispersion, a^, and for 
literature samples of classical radio galaxies and optically selected ellipticals. Only the 
optically selected in-active galaxies are found to exhibit a correlation between M^^ 
and radio luminosity. In contrast, the radio powers of the AGN in the samples do not 
correlate with M^h at all, with objects at a given black hole mass ranging over 7 orders 
of magnitude in radio power. 

We have been able to tie in the local population of powerful radio sources with its 
parent population of in-active elliptical galaxies: the local black hole mass function has 
been determined using the elliptical galaxy luminosity function and the Faber- Jackson 
and Mbh — cr relations. This was combined with the fraction of radio-loud black holes as 
function of M^h^ as determined from the optically selected galaxy sample, to derive the 
local volume-density of radio galaxies and the distribution of their black hole masses. 
These are shown to be consistent with the local radio luminosity function and the 
distribution of black hole masses as found in the radio selected samples, and confirms 
that elliptical galaxies comprise the large majority of the radio-loud population of 
active galaxies. 
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1 INTRODUCTION 

1.1 Galaxy formation and the fundamental plane 

How and when galaxies were formed and how they have 
evolved since, is one of the central issues of contemporary 
astrophysics today. Elliptical galaxies may be the best ob- 
jects to study galaxy evolution, since their optical light is 
not dominated by relatively recent star-formation, unlike 
that of spiral galaxies, and their observational properties 
are a better reflection of their long-timescale cosmological 
development. 



In the 1980s it was shown that the main global prop- 
erties of local elliptical galaxies, velocity dispersion, surface 
brightness and efl^ective radius, are strongly correlated, and 
that they describe a plane in this 3-dimensional parame- 
ter space (Dressier et al., 1987; Djorgovski & Davis, 1987). 
Initially, this 'fundamental plane' (FP) of elliptical galaxies 
drew attention mainly for its use as a distant indicator, to 
measure the Hubble constant and peculiar motions in the 
local universe. However, it was also realised that it implied 
a strong regularity in the process of galaxy formation, in 
particular that the mass-to- light (M/L) ratio varied only 
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slightly over 4 orders of magnitude in galactic mass (Djor- 
govski & Davis, 1987). 

Since the M/L ratio is a strong function of the age of 
the stellar population, the study of the fundamental plane 
as function of rcdsliift is a sensitive diagnostic of the evolu- 
tionary history of galaxies, van Dokkum, Franx and Kelson 
and collaborators, who observed elliptical cluster galaxies 
with the Keck telescope and the HST, located at z=0.39 
(van Dokkum & Franx 1996), z=0.58 (Kelson et al. 1997), 
and z=0.83 (van Dokkum et al. 1998), showed that the M/L 
ratio evolves as A log M/Lb oc —OAz {flm = 0.3,nA = 0), 
putting strong constraints on the formation epoch of clus- 
ter ellipticals, with Zform > 2.8 for a Salpeter initial mass 
function. 

The Mgb or Mg2 indices arc also valuable measures 
of the age of the stellar population, which is often used 
in conjunction with the FP. The Mg2 index is found to 
be correlated with the central velocity dispersion in local 
elliptical galaxies, with only a small scatter (eg. Dressier 
et al., 1987; Djorgovski & Davis, 1987, Jorgensen, Franx 
& Kjaergaard 1996). Ziegler & Bender (1997) showed that 
for cluster galaxies at z=0.37, this correlation is offset by 
<Mgb> ~ — O.4A which can be fully attributed to the lower 
luminosity-weighted ages of these stellar populations at this 
earlier epoch. 

These studies may seem to suggest that, except for per- 
haps some minor details, the evolution of elliptical galaxies 
is well understood. However, various strong biases may still 
be present in the studies described above. The most seri- 
ous bias may be the progenitor bias (van Dokkum & Franx 
1996): present day ellipticals which were spirals at an ear- 
lier epoch, would not have been included in a sample at 
that epoch, biasing the sample to what arc likely to bo the 
oldest galaxies. Furthermore, there is also an environmental 
bias: present studies are focused on cluster galaxies, while 
it is not known how these clusters evolve with cosmological 
epoch. For example, it is far from clear what a local cluster 
such as Coma would have looked like at high rcdshift, and to 
what present-day structures the high-z clusters have evolved 
into. More generally, connecting populations at one redshift 
with their predecessor at high redshift, without making un- 
certain assumptions about the history of star formation is a 
problem. 

1.2 The role of nuclear activity in galcLxy 
formation and evolution 

The place of active galaxies in the general picture of galaxy 
evolution is unclear. Until recently, AGN were seen as rare 
and peculiar objects, which, although interesting for their 
own sake, were disconnected from galaxy evolution issues. 
However, the recent discovery that every galaxy seems to 
have a central massive black hole, the meiss of which is 
closely related to the total mass of (the bulges of) its host 
(Magorrian et al, 1998; Gebhardt et al. 2000), has changed 
this perception dramatically. It is now clear that the for- 
mation and evolution of the central massive black hole is 
closely linked to that of the host galaxy (eg. Silk & Rees 
1998). In addition, it implies that all galaxies may be capa- 
ble of having phases of central activity, meaning that active 
galaxies are just normal galaxies, with their central black 
hole coincidently caught during a period of (high) accre- 



tion. This places active galaxies back in the centre of atten- 
tion of galaxy formation and evolution studies, especially 
since they are relatively easy to select at high redshift. Re- 
cent claims suggest that non-thermal radio emission from 
ellipticals is strongly correlated with the mass of their cen- 
tral black hole, ranging from dwarf ellipticals to powerful 
radio galaxies (Franceschini et al. 1998). This would make 
the connection between quiescent and active galaxies even 
stronger. 

1.3 Radio galaxies as probes of galaxy evolution 

Since active galaxies play such a potentially important role 
in galaxy evolution, it is important to determine their fun- 
damental structural parameters as function of redshift. An 
additional factor is that selection of AGN allows the study 
of galaxy evolution independently of the optical selection 
biases as mentioned above. Smith, Heckman & lUingworth 
(1990) have studied the fundamental plane parameters of 
nearby powerful radio galaxies, and showed that they in- 
deed reside in normal ellipticals. A similar result has been 
found recently by Bettoni et al. (2002). Unfortunately, it is 
almost impossible to use similar objects at higher redshift, 
since their spectral energy distributions have invariably been 
found to be strongly influenced by AGN light, especially in 
the optical (eg. optical/UV alignment effect, McCarthy et 
al. 1987, Chambers et al. 1987). 

This problem cam be avoided by selecting galaxies con- 
taining a very young radio source. In the very early stages 
of radio source evolution, when the radio source is only 100- 
1000 years old and confined to the central hundred parsecs, 
the contribution of the AGN to the global properties of its 
host galaxy is still minimal (Snellen et al. 1996a, b; 1999): 
broadband photometry showed that their optical to near- 
infrared colours are consistent with passively evolving ellip- 
ticals, indicating that their light is dominated by old stellar 
populations. The dispersion in their R and K band Hubble 
diagrams (0.3 mags) is smaller than that for other classes 
of radio galaxies, indicating a homogeneous and well be- 
haved population of hosts (Snellen et al. 1996a,b). In ad- 
dition, their optical spectra typically only show weak emis- 
sion lines but deep stellar absorption lines, again indica- 
tive of them being dominated by old stellar populations 
(Snellen et al. 1999). For example, the archetype young radio 
source, B0108-I-388, is as powerful as Cygnus A. However, 
it has a remarkably low OIII5007 equivalent width of only 
^4A(Lawrence et al. 1996), corresponding to a line lumi- 
nosity about a factor 50 lower than that of Cygnus A. 

That the contamination from young radio-loud AGN 
is only small or even absent may not come as a surprise. 
The young radio source may not yet have pierced through 
the dense and dusty central regions of the host. AGN re- 
lated light is also likely to still be conflned to the inner 
hundreds of light years, and obscured from our view. There 
is indeed much evidence through HI absorption studies of 
young radio-loud AGN that large column densities of gas 
are present (eg. Carilli et al. 1999; Peck, Taylor & Conway 
1999) 

In this paper, we investigate the use of the host galaxies 
of young radio-loud AGN as probes of galaxy evolution. The 
connection between radio activity and black hole mass, as 
infered from their central velocity dispersion, has also been 
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investigated. We present high dispersion spectroscopic and 
optical CCD imaging data of a sample of 7 young radio-loud 
AGN located between 0.06 < z < 0.31. Section 2 describes 
the properties of young radio-loud AGN and the selection of 
the sample. Section 3 describes the observations, reduction 
and analysis. Section 4 and 5 present the results and discus- 
sion. For consistency with previous papers on this subject, 
we use a cosmology with Ho=50 km sec~^ Mpc~^, ftrn=i O, 
^Ia=0, throughout this paper, unless stated otherwise. 



2 YOUNG RADIO-LOUD AGN AND THE 
SELECTION OF THE SAMPLE 

2.1 Characteristics of young radio- loud AGN 

In the course of the last decade, a class of extragalactic ra- 
dio source has been identified as being very young. This 
has opened many unique opportunities for radio source evo- 
lution studies. Unfortunately, the nomenclature and use of 
acronyms in this field of research is rather confusing. This is 
mainly caused by the different ways in which young radio- 
loud AGN are selected. Selections of young sources are 
made in two ways, the first based on their broadband ra- 
dio spectra, and the second based on their compact mor- 
phology. A convex shaped spectrum, peaking at about 1 
GHz distinguishes young radio sources from other classes 
of compact radio sources. In this case they are called Gi- 
gahertz Peaked Spectrum (GPS) radio sources (eg. O'Dea 
et al. 1991, O'Dea 1998). Similar objects, which are typi- 
cally an order of magnitude larger in size, have their spec- 
tral turnovers shifted to the 10 — 100 MHz regime, causing 
thcin to be dominated at cm wavelengths by the optically 
thin parts of their spectra. These are called Compact Steep 
Spectrum (CSS) radio sources to distinguish them from the 
general population of extended steep spectrum sources (eg. 
Fanti et al. 1990). 

On the other hand, young radio sources are found in 
multi-frequency VLBI surveys, in which they can be recog- 
nised by compact jet/lobe-like structures on both sides of 
their central core. They are called Compact Symmetric Ob- 
jects (CSO, Wilkinson et al. 1994). Their double sided 
structures clearly distinguish them from the large majority 
of compact sources showing one-sided core-jet morphologies. 
Larger versions of CSOs are subsequently called Medium 
Symmetric Objects (MSO) and Large Symmetric Objects 
(LSO). 

The overlap between the classes of CSO and GPS 
galaxies is large and we believe that they can be consid- 
ered to be identical objects. However, a substantial frac- 
tion of GPS sources are optically identified with high red- 
shift quasars, which in general show core-jet structures 
(Stanghellini et al. 1997). The relationship between GPS 
quasars and GPS/CSO galaxies is not clear (Snellen et al. 
1998b), and the former may not be young objects et all. 

Although it was always speculated that GPS sources 
were young objects, only recently has strong evidence been 
found to support this hypothesis. Monitoring several proto- 
type GPS/CSO sources over a decade or more using VLBI, 
allowed Owsianik & Conway (1998) and Owsianik, Conway 
& Polatidis (1998) to measure their hotspot advance speeds 
to be ~ 0.1h~^c. These imply dynamical ages of typically 



lO'^"^ years. Similar studies with similar results have now 
also been conducted by Tscliager et al. (2000) and Taylor 
et al. (2000). Additional proof for youth comes from analy- 
sis of the overall radio spectra of the somewhat larger CSS 
sources. Murgia et al (1999) show that their spectra can 
be fitted with synchrotron ageing models, implying ages of 
typically 10^~® years. 

The work of these authors indicate that GPS/CSO 
sources are very young and the likely progenitors of large, 
extended radio sources. However, there may also be a pop- 
ulation of young radio-loud AGN which are just short-lived, 
and do not evolve into large size radio sources (Readhead 
et al. 1994; Alexander 2000). In the latter case, the connec- 
tion of these host galaxies with those of other types of AGN 
would be less clear. 

2.2 Sample selection 

The sample of young radio sources used in this paper is 
based on the sample of bright GPS sources of Stanghellini 
et al. (1998). This complete sample of 33 objects contains 
sources with turnover frequencies between 0.4 and 6 GHz, 
with Ss GHz > 1 Jy, and are located at (5 > —25° and galactic 
latitude |6| > 10°. Of this sample, 19 are optically identified 
with galaxies, of which 6 are found at 2 < 0.31. Two objects, 
B1345+125 and B1404+286, were omitted from the sam- 
ple, since these were known to exhibit broad Seyfert-type 
emission lines in their optical spectra, implying a strong 
contribution of AGN related light. We further excluded 
B0428-f-205 due to its low galactic latitude and high fore- 
ground extinction. 

To enlarge the usable sample and to be able to inves- 
tigate possible correlations between host galaxy properties 
or optical AGN contribution and radio power, we included 
nearby GPS galaxies from the sample of Snellen et al (1995) 
at intermediate flux densities, and from Snellen et al. (1998a) 
at faint flux densities. This involved the only GPS galaxy at 
low redshift from the intermediate sample, B1144-I-3517, and 
two nearby galaxies from the faint GPS sample, B0830+5813 
and B1946+7048. Details of the young radio-loud AGN are 
given in table 1. Column 1 gives the B1950 name, and col- 
umn 2 gives the J2000 position. Column 3 gives the redshift. 
Column 4 gives the R band magnitude, taken from Snellen 
et al. (1996a), O'Dea et al. (1996) and Snellen et al. (1998b). 
These are approximately corrected to the Cousins filter sys- 
tem, and galactic foreground extinction is taken into ac- 
count, which is given in column 5 (Schlegel et al. 1998). Their 
5 GHz radio power, radio peak frequency, peak flux density, 
angular and projected linear size are given in columns 6 to 
10. 



3 OBSERVATIONS, DATA REDUCTION AND 
ANALYSIS 

The observations for this project are divided into three 
parts. 1) High-dispersion spectroscopy, to determine the ve- 
locity dispersions and absorption line strengths of the ob- 
jects. 2) Deep I-band CCD imaging to determine the surface 
profiles of the objects, and 3) B, V and R-band observations 
to determine the rest-frame B-band effective surface densi- 
ties of these galaxies. 
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Figure 1. The part of the spectra used to measure the velocity dispersion, shown in order of increasing redshift. The thick lines shows 
the best fit to the velocity broadened stellar spectra. Note that in some cases, eg. for B0941-0805, the [NIJ5199A emission line had to be 
removed before the fitting procedure. 
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Figure 2. The normalised galaxy profiles with their best fits. The profile of B0830+5813 is found to flatten in the central part. 
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Table 1. Details of the young radio- loud AGN in the sample. 
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Figure 3. Tlic galaxy spectra of which the velocity broadened K-stax spectra are subtracted. The three adjacent panels to the right zoom 
m on the [0111]4959A, [Olll]5007A, and the [NI]5199A emission lines, with the dashed lines indicating their rest-frame wavelengths, 
assuming the absorption line redshifts. 



3.1 High-dispersion Spectroscopy 

The spectroscopic observations wore carried out using the 
4.2-m William Herschel Telescope (WHT) located at the 
Roque de Los Muchachos Observatory, La Palma in Spain. 
Observations were carried out during two observing sessions 
July 29/30, 1997 and November 9-11, 1998. Both observing 
runs were under rion-photorrictric and variable seeing condi- 
tions. Table 2 gives the log of the observations, with in col- 
umn 1 the object name, and column 2 the observing run(s) 
during which the source was observed (l=July '97, 2=Nov 
'98). The exposure times per source ranged from 30 minutes 
to several hours, dependent on the redshift of the target and 
weather ( a significant proportion of the observing time was 
affected by clouds) . 

The red arm of the ISIS long-slit spectrograph was used, 
with the R600R grating, which has 600 grooves/mm. The 



spectra were centred on the redshifted Mg b absorption fea- 
ture near 5170A. The approximate central wavelengths are 
given in column 3 of table 2. The detector used was a Tck- 
tronic 1024 x 1024 CCD, which resulted in a dispersion of 
0.79 A per pixel, and a spatial resolution of 0.3 arcseconds. 
For all observations a slit- width of 1.0 arcseconds was used. 
The spectral resolution for this setup was determined to be 
1.6 A FWHM. The rest-frame spectral resolution for the ob- 
servations of each target is given in column 4 of table 2. The 
resulting signal-to-noise ratio per pixel is given in column 6. 



Reduction and calibration of the spectra 

The data were reduced in IRAF. The individual raw im- 
ages were first bias subtracted. After each target, an ex- 
posure was taken of the Tungsten lamp with the telescope 
pointing in the same position. This was normalised and cor- 
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Table 2. Parameters of the WHT spectroscopic observations 
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rected for illumination effects, and was used to Hatfield the 
science images and to correct them for fringing. The data 
were then background subtracted and combined, and a one- 
dimensional spectrum retrieved. Although the spectroscopic 
observations were not taken under photometric conditions, 
observations of spectroscopy standard stars were used to re- 
move the response of the instrument from the spectra. 

The wavelength calibration was performed using a 
CuAr calibration lamp. The wavelengths of the main sky- 
lines in the non-background subtracted images were checked, 
and if necessary, a small shift in wavelength was applied. 
Finally, the spectra were corrected to the heliocentric stan- 
dard. The uncertainty in the absolute wavelength calibration 
is estimated to be 0.2A (~10 km/sec). 



Mg b absorption line equivalent widths 

The strength of the Mgb absorption feature was measured 
in an 'index passband' flanked by two 'pseudo-continua' (see 
Worthey 1994). The wavelength-ranges of the index and 
continuum band are as defined in Worthey et al. (1994, 
tabic 1). Unfortunately, due to a varying contribution of 
the [NI]5199A emission line in the upper pseudo-contiimum 
band of the Mgb index, the standard Mgb bands could 
not be used. Instead, an upper pseudo-continuum band at 
5205-5220A was used, avoiding this line. Our template 
star spectra were used to calibrate this new measure of 
the Mg b index, and we found the empirical relation Mg b= 
0.677 -I- 0.790xMgbait, to correct these values. 

The measured absorption indices had to be corrected 
for the velocity broadening of the spectra. This effect smears 
out the absorption lines, possibly extending them to outside 
the index passband and into the pseudo-continuum bands, 
resulting in a reduction of the measured absorption index. 
This effect, which is a function of the velocity dispersion, is 
corrected for by estimating its influence by velocity broaden- 
ing the stellar-spectra. The corrections, in the order of 10%, 
were compared to that derived by Trager et al. (1998), and 
found to be similar. 

Since elliptical galaxies have radial gradients in line 
strengths, the effects of different aperture size must also be 
taken into account. Wc used a similar approach to that of 
Ziegler & Bender (1997), using Mg2 profiles analysed by 
Gonzales & Gorgas (1995). The values were corrected to a 
radius of re/4, resulting in an adjustment between —0.25 
and +0.28. 



The velocity dispersions and absorption-line redshifts 

The spectra were analysed by a Fourier Quotient method 
(Simkin 1974; Sargent et al. 1977), using the software pack- 
age IDL. To determine the velocity dispersion and the ab- 
sorption line redshifts of the galaxies, their spectra were 
compared to that of K-stars. These stars proved to provide 
the most appropriate absorption line spectra for compari- 
son to these old elliptical galaxies. For each galaxy, the con- 
tinuum was subtracted using a 5th order polynomial, after 
which the extremities of the spectrum were tapered using 
a Cosine-bell function. The Fourier transform of this spec- 
trum was divided by that of the continuum subtracted and 
tapered star spectrum. The result represents the broaden- 
ing function which was fitted by a Gaussian, and gives the 
redshift, velocity dispersion, and relative strength of the ab- 
sorption lines of the galaxy to that of the template star. 
This fitting in Fourier space was typically performed over a 
wave number range corresponding to 3-lOOA. Variations in 
redshift and velocity dispersion between different template 
stars were found to be minimal. In some cases the procedure 
was preceded by the removal of a small contribution from 
the [NI]5199A emission line. 

Since elliptical galax;ies exhibit radial gradients in their 
velocity dispersions, the measured velocity dispersion had 
to be corrected for the effective aperture size of the obser- 
vations. Wc used the relation found by Jorgcnscn, Franx & 
Kjaergaard (1995) to correct the velocity dispersions to an 
aperture diameter of 1/4 of the effective radius, re/4. This 
resulted in aperture corrections up to 20 km sec~^. 



The emission line equivalent widths and redshifts 

The best fitting velocity broadened stellar spectrum was 
subtracted from each galaxy to isolate the emission lines 
and measure their equivalent widths. The velocity shifts of 
the emission lines were determined by comparing the mear- 
sured peak wavelengths with that of laboratory rest-frame 
measurements. 



3.2 I-band CCD imaging 

The imaging observations to measure the effective radii of 
the galaxies were also carried out using the WHT, this dur- 
ing the second observing session at November 9-11, 1998. 
The weather was non-photometric and the seeing condi- 
tions varied considerably. The auxiliary port camera at the 
Cassegrain focus was used with a Tektronics 1024x1024 
CCD detector, which resulted in a 0.11" pixclsizc and an 
unvignetted (circular) field with a diameter of 1.8'. After ev- 
ery exposure of a few minutes, the telescope was offsetted by 
tens of arcseconds, to allow more accurate sky-subtraction. 
Total exposure times ranged from 10 to 45 minutes per 
galaxy. 

The data were reduced in IRAF. The individual raw im- 
ages were first bias subtracted and flat flelded. These were 
then combined using median filtering to produce a fringe- 
images which was then appropriately normalised and sub- 
tracted from the frames. The centroids of bright stars on the 
frames were used to determine their relative shifts, which 
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Table 3. The results of the imaging and spectroscopic observations. 
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11.9±1.7 


28.0±0.8 


25.11 


4.33 


3 


-110 


B0941-0805 


0.22790±0.00012 


148±17 


2.9±0.3 


14.0±0.9 


24.03 


2.73 


88 


-1-20 


B1144+3517 


0.06313±0.00013 


202±18 


4.9±0.5 


8.3±0.5 


21.82 


4.11 


25 


-5 


B1819+6707 


0.22071±0.00011 


221±16 


4.0±0.5 


18.8±1.5 


22.85 


4.29 


6 


+100 


B1946+7048 


0.10083±0. 00009 


171±20 


7.0±0.3 


17.7±0.5 


24.33 


4.64 


8 


-0 


B2352+4933 


0.23790±0.00016 


201±17 


1.6±0.3 


7.9±1.0 


22.15 


4.35 


35 


-1-350 



were then combined to produce the final image. No flux cal- 
ibration was attempted due to the non-photometric condi- 
tions. 

The effective radii 

The reduced I band images were used to the determine the 
eflfective radii of the galaxies. This analysis was done in the 
IDL software package. Firstly, the centroid of a galaxy was 
determined and its ellipticity fitted at an isophote at 20% 
of the value of the central peak. The level of intensity was 
measured for an increasing major axis, using the whole area 
of the ellipse, in this way making full use of the data. Some 
areas of the image, including nearby objects, were excluded 
from the analysis. The central parts of the galaxies, which 
were influenced by the variable seeing, were also excluded. 
The resulting profile was fitted by a r^''* de Vaucouleur pro- 
file. 

3.3 B,V and R-band CCD imaging 

Short CCD exposures of the objects were taken with the 
Jacobus Kapteyn Telescope, which is also located at the the 
Roque de Los Muchachos Observatory, La Palma in Spain. 
These served to determine their surface brightnesses at their 
effective radii in rest-frame B. Observations were carried 
out during an observing session from June 14 to June 23 
2001, and during service observations on July 4 and August 
19 2001. All JKT observations were under good photomet- 
ric conditions, and typical seeiing of 1 — 1.5". We used the 
CCD camera with the default SITe2 detector, which has 
2048x2048 24/im pixels, resulting in an image scale of 0.33 
arcsec/pixel. We observed the objects in B and V , or in V 
and R filter (depending on the redshift) for 2x300 sec in 
each filter. Standard star fields from Jorgensen (1994) were 
observed throughout the nights at several zenith distances. 

The data were reduced in a standard way, in a similar 
way as the WHT I-band data. The photometric zero-points 
were found to vary by less then < 0.1 magnitude between 
the different observing nights. 

Rest-frame B-band effective surface brightness 

To determine the effective surface brightness, le in rest- 
frame B band, first the effective surface brightnesses were 
measured in each filter: the ellipticity and position angle of 
the galaxies as determined from the deep I-band observa- 
tions were used to determine the surface brightness as func- 
tion of radius for each galaxy, excluding some areas contain- 



ing neighbouring objects. This ID profile was then fitted 
with an r^^** law, using the effective radii as found from the 
I-band observations, allowing 7e and the sky-level to vary. 
As a consistency check, the total intensities of the galaxies 
were also measured, using a large aperture of >> re- For 
a perfect 2-dimensional De Vaucouleur profile, these relate 
as hot = 22.66/ere. For the galaxies free of nearby objects, 
they were found to match within 0.1 magnitude. 

The V and R band fiux density scales were corrected to 
match the B band zero-points, adopting an absolute spec- 
tral irradiance for a zero magnitude object of 6.61 x 10~^, 
3.64 X 10"^ and 2.25 x 10"'' erg cm"^ sec"^ A"^ in the 
B, V and R band respectively, and were also corrected for 
galactic foreground extinction using the values of Schlcgel et 
al. (1998). Interpolation between the surface brightnesses in 
the different filters was used to obtain the rest-frame B-band 
surface brightness. This was then corrected for the cosmo- 
logical dimming and bandpass effects of (1 -|- zf . 

For two galaxies, B0830+5813, and B0914-0805, only 
R band images were available. For these we applied a K- 
correction, derived from the broadband colours of the other 
galaxies in the sample. 



4 RESULTS AND DISCUSSION 

The results of the imaging and spectroscopy are given in 
table 3. Column 1 gives the B1950 name, columns 2 and 
3 give the redshift and the velocity dispersion corrected to 
an aperture of re/4. Column 4 and 5 give the effective radii 
in arcseconds and kpc. Column 6 gives the brightness level 
at the effective radius for rest-frame B-band, corrected for 
cosmological dimming and galactic foreground extinction. 
Columns 7 and 8 give the rest-frame equivalent widths for 
the Mg b band (also corrected to an aperture of re/4) and 
the [OIII]5007A emission line. Column 9 gives the veloc- 
ity shift of the [OIII] line with respect to absorption line 
redshift, where a positive value indicates a redshift. The un- 
certainties in Je are typically 0.05 magnitude, while the un- 
certainties in the equivalent widths are ~ 10%. 

The part of the spectra used to meeisure the velocity 
dispersions of the galaxies is shown in figure 1, in order of in- 
creasing redshift. The thick lines indicate the best fit to the 
velocity broadened K-star spectra. The measured velocity 
dispersions range from from 146 km sec~^ for B0830-)-5813, 
to 228 km sec"^ for B0019-000. Figure 2 show the nor- 
malised galaxy profiles with their best fits. The observed 
effective radn range from 1.1" for B009-0001, to 11.9" for 
B0830-I-5813. Note that the profile of B0830-I-5813 sUghtly 
flattens in its inner part. The best fltted velocity broadened 
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Figure 4. The fundamental plane of the host galaxies of young radio-loud AGN (solid squares), with on the horizontal axis their effective 
radii, and on the vertical axis a combination of the velocity dispersion and their effective brightness, corrected for the cosmological 

evolution in the mass-to-liglit ratio as found by van Dokkum ct al. (1998). The solid line indicates the fundamental plane of elliptical 
galaxies, as derived from nearby galaxies by Faber et al. (1989; small circles). The diamonds are the radio galaxies measured by Smith 
et al. (1990). 



K-star spectra were subtracted from those of the galaxies 
to obtain their emission line spectra. These are shown in 

figure 3. The throe adjacent panels to the right zoom in 
on the [OIII]4959A, [OIII]5007A, and the [NI]5199A emis- 
sion lines, with the dashed lines indicating their rest-frame 
wavelengths, assuming the absorption line redshifts. it shows 
that the emission line redshifts correspond well with those 
infered from the absorption lines, except for B2352+4933, 
which emission lines axe clearly redshifted by 350 km/sec. 



4.1 The Fundamental Plane 

The main aim of the project was to investigate whether 
young radio-loud AGN could be used to study the cosmolog- 
ical evolution of the fundamental plan for ellipticals, being 
selected independently of optical selection biases. Figure 4 
shows the fundamental plane with the angle of projection 
in the a, le, re - space chosen to be similar to that used in 
previous studies, with on the horizontal axis the logarithm 
of the effective radius, re (in kpc), and on the vertical axis 

1.2 times the logarithm of the velocity dispersion a (in km 
sec~^) minus 0.83 times the logarithm of the brightness level 
at the effective radius, le, with le = 10"° ''''^. The data are 



0.6 I 

J 0.4 . 




c 

3 - ■ 

-0.4 7 . 

o 

-0.6 r I I I ' 

0.0 0.1 0.2 0.3 0.4 

Redshift 

Figure 5. The offset of the host galaxies to the local fundamental 
plane for cluster ellipticals as function of redshift. The dotted line 
is the evolution in mass-to-light ratio as found by van Dokkum 
et al. (1998). The symbols are as in figure 4. There is a possible 
hint that the more distant galaxies are brighter, but more data 
are needed to confirm this. 
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Figure 6. The Mgb-a- relation for the host galaxies of young 
radio- loud AGN (squares). The solid line indicates the Mgb-cr 
relation for local ellipticals from Faber et al. (1989; small circles). 



0.5 



b -0.5 - 



-1.0 - 



-1.5 





















* 



20 40 60 80 100 

OIIIj^QQ.-, equivalent width [A] 

Figure 7. The offsets of the Mg b indices to the Mg b-cr relation 
for local ellipticals by Faber et al. (1989), as function of [OIII] 
equivalent width. The correlation indicates that the Mgb index 
has likely been influenced by continuum contamination from AGN 
related light in the nucleaj: spectra. 



corrected for the cosmological evolution in mass-to-light ra- 
tio as found by van Dokkum et al. (1998), corresponding 
to Amag = 1.22X2;. The solid squares show the galaxies in 
the sample, and the diamonds axe the radio galaxies from 
Smith, Heckman & Illingworth (1990). The small circles are 
the local ellipticals from Faber et al. (1989) with heliocentric 
velocities >1500 km/sec, with the solid line their best linear 
fit. The galaxies in our sample follow the fundamental plane 
well, with an average offset of 0.01 and a dispersion of 0.17, 
similar to that of the Smith et al. sample, but with a slightly 
higher dispersion than the local galaxies in the Faber et al. 
sample (0.13). 

To investigate whether evolution in M/L is visible over 
the small redshift range in our sample, the offsets of the dat- 
apoints to the local fundamental plane, AFP, as function of 
redshift is shown in figure 5. The symbols arc as in figure 4, 
with the grey band indicating the dispersion in the local fun- 
damental plane, and the dashed line showing the evolution 
as found by van Dokkum et al. (1998) for cluster ellipti- 



cals. The galaxies containing young radio sources located at 
z>0.2 seem to have on average a slightly positive offset from 
the local fundamental plane, with AFP=0.10±0.06. There- 
fore evolution is only detected at a 1.7ff level in this small 
datasct. Clearly, more data arc needed, preferably out to 
higher redshift, to confirm evolution, and to test the find- 
ings by van Dokkum et al. 



4.2 The Mgb - a relation 

The absorption line strengths of a galaxy also give a good, 
independent insight in its stellar population. The Mgb-cr 
relation, and its (Icpoudcnce on redshift, has been shown 
to be a powerful tool to study galaxy evolution (Ziegler & 
Bender 1997). The Mgb-cr relation for the galaxies in our 
sample is shown in figure 6. The small circles arc the data 
from the local Faber et al. sample, where the best fit is indi- 
cated by the solid Une. These datapoints were converted Mg2 
magnitudes, using the relation MgbA~^=15.0Mg2 (Ziegler 
& Bender 1997). The Mgb index of B0941-0805 is clearly 
lower than expected from the local relation. Since this ob- 
ject exhibit a strong emission line spectrum, we suspect that 
this Mg b deficit is caused by a contamination of the contin- 
uum by AGN related light. In figure 7, the offsets of the Mg b 
indices to the local Mgb-cr relation is shown to be clearly 
correlated with OIIIsoot equivalent width, indicating that 
AGN contamination indeed plays a role. The Mg b deficit in 
B0941-0805 can be explained by an AGN continuum con- 
tribution of 35%. For the other galaxies in the sample the 
AGN contribution is less than <5-10%. Note that this AGN 
contamination is for nuclear spectra, using a narrow slit. It 
is unlikely to have any effect on the determinations of the 
surface brightness profiles as done at and beyond re, and 
it is therefore improbable to affect the determination of the 
fundamental plane parameters. However, this AGN contri- 
bution, although generally small, makes investigating galaxy 
evolution using the Mg h-a relation not very useful for radio 
galaxies. 



4.3 Black hole masses and their connection with 
radio power 

The work by Gebhardt et al. (2000) and Ferrarese & Merrit 
(2000) have shown that the mass of a central black hole is 

strongly correlated with the central velocity dispersion of 
its host galaxy. According to Gebhardt et al. the relation 
follows 



Mbh = 1.2 X 10" X 



200 km/s 



which is found to have a lower dispersion than the correla- 
tion with absolute (bulge) magnitude. This Mth — cr relation 
allows us to infer black hole masses from the central veloc- 
ity dispersion measurements, for the galaxies in our sample, 
and those in Smith et al. and Faber et al. Uncertainties in 
the black holes masses are dominated by the scatter in the 
Mbh — o relation, which is estimated to be a factor of ~2 
(Gebhardt et al. 2000). 

A variety of authors have discussed the relation between 
Mbh and radio power. Long before the correlations between 
black hole mass and central velocity dispersion or galactic 
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Figure 8. The radio luminosity at 1.4 GHz versus black hole mass, as derived from the central velocity dispersion. The squares axe those 
objects selected in the radio (from this paper and Smith et al. 1990). The open and closed circles are optically selected galaxies from the 

Fabcr ct al. (1989) sample, with the former being point sources in the NVSS, and assumed to be radio-quiet, and the the latter being 
extended in the NVSS and assumed to be radio-loud AGN. The few diamonds at low black hole mass are galaxies added from the local 
group. The solid line is the Mi,;j — Pradio relation as found by Pranceschini et al. (1998). 



bulge mass were established, evidence has been provided of 
a relation between galactic bulge luminosity and/or velocity 
dispersion with radio luminosity L_r (Hcckman 1983, Nelson 
& Whittle 1996). Pranceschini et al. (1998) showed that for 
a small sample of nearby galaxies, the non-thermal radio 
emission is strongly correlated with the black hole mass, 
from the dwarf elliptical M32 to the powerful radio galaxy 
M87, with a surprisingly small scatter, which was confirmed 
by McLure et al. (1999). More recently, other studies (eg. 
Laor 2000; Ho 2002) show that AGN are in general located 
above the overall Mbh^LR relation as found by Franccschini 
et al. Interestingly, Dunlop and McLure and collaborators 
(Dunlop et al. 2002; Dunlop & McLure 2002) have suggested 
that there is an Ma^-dependent lower and upper limit to the 
possible radio output of a black hole, ofFsetted from each 
other over several orders of magnitude. 

To investigate these issues further, and in particular the 
relative position of radio galaxies on the Adth — Lr relation, 
we correlated the galaxies from Faber et al. with the NVSS- 
VLA 1.4 GHz survey (Condon et al. 1998) and determined 
their radio luminosities and sizes. 

The Mbh — Lr relation for the galaxies in the Faber et 
al. (circles) , and the radio selected galaxies from Smith et al. 



and this paper (squares) axe shown in figure 8. The galaxies 
in Faber et al. which exhibit extended radio structure, and 

therefore clearly identified as active galaxies, arc indicated 
by solid circles. A few objects from the local group are added, 
and are indicated by diamonds. 

To avoid the diagram being crowded by meaningless 
upperlimits, only those objects are plotted, if their radio 
flux would have been detectable in the NVSS if they ra- 
diate at a luminosity expected from the Franceschini rela- 
tion. Although this produces a varying redshift cut-off with 
black hole mass, it does not bias the result towards par- 
ticular radio luminosities. Clearly, the radio luminosities of 
the presumably inactive galaxies correlate well with black 
hole mass, although they seem to be slightly fainter than 
expected from the Franccschini relation, and have a larger 
scatter, with the radio power of inactive galaxies at a given 
black hole mass ranging more than 2 orders of magnitude. 
The active galaxies (those selected in the radio, and those 
selected from the Faber et al. sample exhibiting extended 
NVSS emission), do not show any sign of a correlation with 
black hole mass. At a given Mih the radio luminosity can 
span over 7 orders of magnitude, making radio power un- 
suitable as a black hole mass estimator. We note that the 
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Figure 9. The fraction of radio-loud AGN as function of black 

hole mass, as derived from the optically selected sample of 
Faber et al. (1989). The distributions are explained in the text. 
The smooth curve is the function P(AGN — Mth) used to deter- 
mine the density and M;,^ distribution for radio-loud AGN with 
f'l.4GHz>23.0 W/Hz. 
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Figure 10. The distribution of black hole masses for the radio 

selected galaxies. The smooth curve (arbitrarily scaled) shows 
the predicted distribution of Mi,h derived from the local elliptical 
galaxy luminosity function and the fraction of black holes with 
powerful radio sources as function of Mj^i as derived from the 
optically selected sample. 



black hole masses for radio galaxies as derived by Dunlop 
et al. (2002) are on average about a factor of 4 higher than 
given hero (assuming Mi,;,=0.0013Mbuige). This is caused by 
their different method used to determine Mbh, based on the 
bulge luminosity — bulge mass — black hole mass relation: 
the average absolute R-band magnitude for the Smith et 
al. sample (—23.71 ± 0.12, assuming V-R=0.6) is the same 
as that for the radio galaxies in the Dunlop et aJ. sample 
(—23.66 ± 0.16), indicating that the properties of the host 
galaxies in the two samples are the same. Note that the 
differences between these methods result mainly in an un- 
certainty in the absolute values of Mi,h , not in their relative 
values. All the radio galaxies are within the allowed range 
of radio power for their black hole masses as proposed by 
Dunlop et al. (2002). 

4.4 Connecting the local populations of radio 
galaxies and inactive ellipticals 

Since the radio galaxies are selected virtually independently 
of their optical properties, and the optically selected galax- 
ies independently of their radio properties, the samples can 
be used to connect the two populations. Figure 8 shows that 
there is a clear difference between the Mbh distribution for 
radio selected and optically selected AGN, with those se- 
lected in the radio, being biased towards lower black hole 
masses. We believe this is the natural consequence of the 
Faber et al. sample strongly favouring intrinsically more lu- 
minous galaxies, caused by its apparent magnitude limit. In 
contrast, the selection of objects on radio emission is practi- 
cally volume based, with the parent population of potential 
radio-loud AGN containing relatively many more galaxies 
with lower velocity dispersions and black hole masses than 
the Faber et al. sample, which causes the sample of radio- 
selected AGNs to cover more galaxies with low Mth- 

These black hole mass distributions form the key to un- 
derstand the connection between the populations of radio- 
loud and radio-quiet galaxies. First we determined the frac- 



tion of radio-loud AGN in the Faber et al. sample as a func- 
tion of black hole mass. This is shown in figure 9, with the 
solid and dotted lines indicating the fraction of galaxies with 
PiAGHz > 23.0 and > 24.0 W/Hz respectively, the former 
decreasing from about 1/5 at Mbh = 10^ M0 to < 1/20 at 
Mbh = 10^"''. Note that a different subsamplc has been used 
for this than in figure 8, with such a distance limit that 
it can be determined from the NVSS whether the galaxy 
has a radio power above or below the cut-off. The dashed 
line shows the fraction of objects as function of Mbh with 
a radio power more than 100 times that expected for the 
Mbh — f radio relation by Franceschini et al (1998). It shows 
that although the radio selected AGN do not have black 
holes with M-bh < 3 x lO^'Afg, black holes below this mass 
can still be active but only with a very low probability. Sev- 
eral objects with Mbh ~ IO'^Mq have been found with ex- 
tended double-lobed low power radio sources in the Faber 
et al. sample. 

With this estimate of the fraction of black holes produc- 
ing powerful radio sources as a function of Mbh, the popu- 
lations of radio-loud and radio-quiet galaxies can be linked. 
The absolute magnitude of a galaxy is connected to its cen- 
tral velocity dispersion by the Faber-Jackson relation (Faber 
& Jackson 1976), which in itself is related to the black hole 
mass via the Mbh — o" relation. In this way, the local dis- 
tribution of black hole masses can be roughly deduced from 
the local galaxy luminosity function. The combination of the 
Faber-Jackson relation, which was re-determined for Ue from 
the galaxies in the Faber et al. sample, and the Mbh — fe 
relation, results in a local black hole mass distribution of 

A* / i\/r \°-'^'^ /a/t \0-65(i+«) 

= o..54exp(-M^) dM„ 

with Log M* = -(8.818 + MS)/1.616, where M%, a, and 
(j)* arc the Schcchter Function parameters of the local lumi- 
nosity function in B band. Our estimate of the fraction of 
black holes producing powerful radio sources should now 
link this distribution of Mbh to the local radio luminos- 
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ity function. For our calculation we assume that the frac- 
tion of radio-loud black holes with Pi.4ghz > 23.0 W/Hz is 
0.25 X (1 -I- {Mbh/1-5 X 10® Mq)^), which is represented by 
the smooth curve in figure 9. Using the Schechter param- 
eters for the local luminosity function of elliptical galaxies 
only, as derived from the 2dF Galaxy Redshift Survey (type 
1 gals; Madgwick et al., 2002) we derive a local radio source 
density of 1.6 x 10"^ Mpc"^ at Log Pi.4GHz > 23.0 W/Hz. 
This corresponds very well to the density of radio sources 
above this luminosity cutoff as derived from the local radio 
luminosity function, which is 1.3 x 10~^ Mpc~^ (eg. Dunlop 
& Peacock, 1990). In addition, the distribution of the black 
hole masses in the radio selected samples is shown in Figure 
10. Although these samples arc not complete, we believe that 
they are roughly selected at random from the general popu- 
lation of radio sources with Log Pi.4GHz > 24.0 W/Hz. The 
smooth curve represents the expected distribution of black 
hole masses from the calculations above, which is strikingly 
similar. These results confirm that elliptical galaxies must 
comprise the large majority of the radio-loud population of 
active galaxies. 



5 SUMMARY 

We have determined the fundamental galaxy parameters for 
a sample of 7 host galaxies of young radio-loud AGN at 
z<0.3, and have shown that they follow the fundamental 
plane of elliptical galaxies. There is some indication that 
the more distant objects in the sample have higher M/L 
ratios, as expected from stellar evolution. As these galaxies 
with young radio sources are relatively easy to select, and 
the contributions of the AGN light to the optical spectra are 
minimal, they can readily be used to study the evolution of 
the fundamental plane out to z=l, independent of optical 
selection effects. 

Black hole masses are determined for the young radio- 
loud AGN, and for samples of classical radio galaxies and 
optically selected ellipticals, using the Mbh — cr relation. We 
confirm the relation between radio power and black hole 
mass as found by Franceschini et al. (1998), but only for rel- 
atively passive ellipticals. Active galaxies do not show such 
a correlation and can range over 7 orders of magnitude in 
radio power at a given Mt,h- 

The local galaxy luminosity function was combined with 
the Fabcr-Jackson relation and the Mbh — o" relation to ob- 
tain the local distribution of black hole masses. This, com- 
bined with the fraction of black holes with powerful radio 
sources as function of Mi,h, which was determined from the 
optically selected sample, results in a radio source density 
and radio-loud black hole mass distribution. These arc found 
to be consistent with the local radio luminosity function and 
the distribution of black hole masses in the radio-selected 
samples, and confirms that elliptical galaxies must comprise 
the large majority of the radio-loud population of active 
galaxies. 
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